Deterministic Optical Fock State Generation 
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We present a scheme for the deterministic generation of iV-photon Fock states from N three-level 
àtoms in a high-finesse optical cavity. The method applies an external làser pulse that generates 
an 7V-photon output state while adiabatically keeping the atom-cavity system within a subspace 
of optically dark states. We present analytical estimates of the error due to amplitude leakage 
from these dark states for general N, and compare it with explícit results of numerical simulations 
for N < 5. The method is shown to provide a robust source of iV-photon states under a variety 
of experimental conditions and is suitable for experimental implementation using a cloud of cold 
àtoms magnetically trapped in a cavity. The resulting iV-photon states have potential applications 
in fundamental studies of non-classical states and in quantum information processing. 



I. INTRODUCTION 



The generation of non-classical states of light has been central to the confirmation and elucidation of the quantum 
theory of radiation. Today, such work takes on an added importance as part of the advancing field of deliberate 
quantum state engineering, motivated in part by applications in quantum communication and computation. For 
instance, deterministically generated single-mode single-photon states could greatly advance the efficiency and security 
of quantum cryptography. Q They are also a crucial resource for implementing a quantum computing scheme using 
linear òptics, single-photon states, and photodetection g|. These immediate potential applications have spurred 
t^J- \ the development of devices which can produce single photons on demand, such as sòlid state devices which use the 
çT) • Coulomb interaction between strongly confined electrons to produce single exciton states which then decay optically 
j|, ||, q], or devices in which the fluorescence from single, isolated and optically-pumped molècules is collected @, |[. 
The stream of pulses produced from each of these devices has been shown to contain either zero or one, but rarely 
more than one, photon per pulse, thus differing radically from a classical coherent-state pulse which would contain a 
Poisson distribution of photon number. Furthermore, while these recent devices all produce fluorescence from a single 
optical cmitter on-demand, none outputs this fluorescence into a practical single mode of the optical field. 

It has been shown recently that the effects of cavity quantum electrodynamics ÍCQED) can in principle be used 
^p" 1 ' to overcome this limitation and produce single-mode single photons on demand || [UJ. In such a scheme, a single 
three-level atom is induced to fluoresce with high probability into a single resonant mode of a high-fincssc optical 
cavity. Experimental evidence for such cavity- induced Raman transitions has been obtained [ p"ï| , and a variant of this 
scheme has very recently been used to generate a sequence of single photons Q. By creating a highly controlled, 
q-i single-quantum-level interaction of àtoms and light, CQED can be used quite generally to produce non-classical states 
JL^ ■ of the electromagnètic field in a single cavity mode. 

In this paper, we analyze a particular desired function of a CQED device that goes beyond the production of single- 
photon states, namely, the deterministic production of a Fock state of the electromagnètic field containing an exact 
number of photons (N) . Such non-classical states are of interest for fundamental tests and applications of the thcories 
of quantum òptics (such as quantum state tomography, as performed recently on the single-photon state |Ï2|| ), as a 
resource for Heisenberg-limited quantum measurements made possible by the production of two orthogonal iV-photon 
states |Ï4| [Ï5[ , and quite generally as a starting point for the controlled engineering of more complex quantum states. 
Starting with a precisely counted number of N àtoms trapped in the confines of a high-finesse optical cavity, we 
consider a scheme in which a classical pump field is ramped up to induce deterministic Raman emission into a single 
cavity mode by each of the trapped àtoms, resulting in an optical field of exactly N photons that are emitted from the 
cavity in a single pulse. The atom-cavity system is constrained to remain within a subspace of optically dark iV-atom 
states, resulting in a high fidelity of production. The present scheme provides a generalization of that proposed for 
the production of single photons from single àtoms || [ÏÏ], [Ï2| and indicates a systematic route to generation of other 
non-classical states. 

Several other theoretical and experimental works have discussed the use of high-finesse cavities for the quantum 
engineering of mesoscopic non-classical optical states. The possibility of producing both Fock states and arbitrary 
coherent superpositions of these inside a cavity by exploiting adiabatic transfer of atòmic ground-state Zeeman co- 
herences in single àtoms was already explored in Ref. A scheme for producing an arbitrary quantum optical 
state using a single two-level atom in a high-finesse cavity has also been presented by Law and Ebcrly, based on the 
arbitrary real-time control of a classical pump field and the coupling to a cavity field |Ï6[ . Experimental evidence for 
Fock states of a microwave cavity field has recently been obtained as a dynamical equilibrium for a stream of Rydberg 
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àtoms passing through a micromaser |Ï7fl . Another approach toward the construction of a Fock state was proposed 
in Ref. ifïst in which a Rydberg atom with a Stark-tunable level splitting is used to transfer photons one-by-one from 
a classically-populated cavity field to another initially empty cavity field. Initial experimental steps towards this goal 
have been recently demonstrated [|Ï9| . 

In contrast to the schemes of Refs. 16, ï?], [ï^] which require delicate temporal control of the atom-cavity coupling, 



our present schcmc yiclds the dcsircd A-photon state for quite arbitrary temporal profiles of the classical pump field. 
This is achieved through a ràpid adiabatic passage which transfers the initial "non-classical" state of N àtoms to 
the non-classical state of A photons within a short period of time. Thus, the quantum nature of the photon field is 
already guaranteed by the initial state of the àtoms trapped in an initially empty cavity. The initial conditions can 
be achieved by lowering a cloud of cold àtoms into the cavity, and pulse control is then entirely contained in the time 
dcpcndcnce of the ramping field. This avoids the need to make use of atòmic motion in controlling the coupling to the 
cavity field, as was required in Ref. || , and is one of the key elements allowing deterministic production of A-photon 
states to be achieved. The number of àtoms in the cavity can be determined by detection of excited state àtoms 
at the single-atom level, as a result of the cavity-mediated shifts of atòmic levels. Such single-atom detection has 
already been demonstrated f2Cfl , and can be readily extended to larger numbers of àtoms. Consequently, the present 
scheme opens the way to deterministic generation of more complicated quantum states of light by first producing 
non-classical states of trapped àtoms (such as spin-squeezed states produced through interatomic interactions [f2~ï] |22j 1 
or by measurement [^3]), and then transferring that state onto the optical field using CQED. 

The remainder of the paper is constructed as follows. A brief review of the deterministic singlc-photon generation 
schemes of Refs. J9|, |Ï0| is given in Section 0, which establishes some common concepts with the present work. A 



discussion of the A-atom/cavity system in a single-mode external field follows in Section III. We demonstrate here 



the existence of a family of optically-dark coupled A-atom/cavity states and show generally how adiabatic ramping of 
an external field can be used to generate A-photon emission from the cavity. Detailed analysis of the energy spectrum 
of the closed A-atom/cavity system as a function of the ramp time provides estimates of the populations in the 
bright states and also leads to analytic estimates of the energy gap required for limiting adiabatic state transfers. In 
Section |y| we then present analysis of the cavity decay responsible for the A-photon emission, treating in detail the 
cffccts of spontancous emission and non-adiabaticity on the output states. We obtain analytic estimates of error rates 
deriving from these contributions that scale linearly in the total number of àtoms, resulting in a constant relative error 
in the output A-photon states and guaranteeing production of an A-photon state with high fidelity. In Section [v| we 
present numerical simulations of the open system using the quantum jump approach. These numerical calculations 
are used to explore the sensitivity of the scheme to critical experimental parameters, as well as to explore the limits of 
our analytical estimates of the error bounds. Finally, in Section VI we summarize and indicate directions for further 
work and for experimental implcmcntation. 



II. ATÒMIC STATES AND DETERMINISTIC SINGLE PHOTON GENERATION 



Our A-photon generation scheme uses the same infernal atòmic structure used in the single-photon proposals of 
Refs. || [To| , namely àtoms having three infernal levels, labeled |0), |1) and |2), arranged in a A configuration such 
that states |0) and |2) are non-radiating atòmic states while |1) is an excited state connected to states |0) and |2) 
by allowed transitions (FigureJ^). Levels |0) and |1) are typically hyperfine levels of the electrònic ground state. We 
assume that the allowed transitions to state |1) can be addressed selectively. This can be achieved, e.g., as a result 
of polarization selection rules, or merely due to a large energy difference between states |0) and |2). The three- level 
àtoms are located in a high-finesse optical cavity which supports a resonant mode having vacuum Rabi frequeney g 
(i.e., the Rabi frequeney due to the presence of a single photon in the cavity mode) that couples states |1) and |2). 
The cavity mode has frequeney lu c which can in general be detuned from the atom resonance W12 by A = w c - lüi2 
(see Figure ||) . 

In the schemes of Refs. Q and [ fÏ0| , single-photon generation is accomplished by exposing a single atom in infernal 
state |0) to a classical làser field of frequeney uj r and Rabi frequeney r(í), which is controlled dynamically. The làser 
frequeney, u r , is chosen to be resonant with the cavity-mediated Raman transition between states |0) and |2), i.e., 
Lú r — lü c = lü2q or u) r — ujiq = A. This làser conneets the states |0) and |2) through a cavity-mediated Raman transition 
and induces the fluorescence of a cavity photon by the atom. The photon èxits the cavity into a single cavity-output 
mode, and thereby generates the desired single-photon state. Refs. |j| and jïo| showed that under suitable conditions 
on the external pulse field, the single photon can be emitted deterministically. 

To understand the operation of such a deterministic device, it is helpful to first consider the atom/cavity system as 
a closed quantum system, i.e., we ignore the decay of cavity photons to cavity-output modes that actually produces 
the desired Fock state outside the cavity, as well as the possible spontaneous decay from the excited state |1) to modes 
other than the cavity mode. Spontaneous emission lowers the fidelity of deterministic photon generation, and clcarly 
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needs to be avoided or at least minimized. We write the basis states of the single atom/cavity system as |0,0), |1,0) 
and |2, 1), where the first index refers to the atòmic state and the second index gives the number of photons in the 
cavity, i.e., the cavity field is implicitly assumed to be quantized. The interaction Hamiltonian of this closed system 
is given by 



H 



o — 



r 
r — A g 
g 



(1) 



where r(í) ~ no\E(t)/2 is the time-dependent coupling to the external (classical) làser pump field, E(t). Unless 
essential for the analysis, we shall omit the explícit time dependence of r to streamline the notation. As discussed by 
Kuhn et al. ]Ï0| , the dynamics of this system are governed by the presence of a null-valued eigenstate l^o), which is a 
"dark state" containing no population in the excited atòmic state |1), and which is therefore immune to spontaneous 
decay. This dark state exists for all vàlues of r(t) and is given by 



l*o> = 



1 



9 Z 



(9\0)-r\2)). 



(2) 



It is the presence of this dark state that allows the high-fidelity generation of a single photon in response to a suitable 
choice of r(t). The single-atom/cavity system starts initially in the state 1 0, 0) , which is the dark state for the initial 
condition r(0) = 0. During a sufficiently slow ramp of r(t) (i.e., a broad pulse of the classical pump làser), the 
atom-cavity system can adiabatically follow the dark state, Eq. (g). For sufficiently large vàlues of r(t) (r 3> g), 
|\&o) ~ |2, 1), i.e., a single cavity photon is produced with a high degree of certainty. This cavity photon then rapidly 
decays from the cavity, resulting in a deterministic single-photon source that acts within a time interval specified by 
the period of the external pump. Law and Kimble (^] and Kuhn et al. JÏ0| have presented numerical calculations to 
assess deviations from this ideal behaviour caused by spontaneous emission, and have explored the extent to which the 
photon emission probability can be controlled by modifying the trigger pulse. Recent experiments by Kuhn, Hcnnrich, 
and Rempe have produced single photons by this method Q. 



III. DETERMINISTIC 7V-PHOTON GENERATION: CLOSED SYSTEM ANALYSIS 

For the generation of arbitrary Fock states of the electromagnètic field, i.e., with arbitrary large photon number 
N, we now consider N such three-level àtoms confined within the optical cavity. We assume that the àtoms are 
indistinguishable in all respects. This provides an important experimental simplification relative to other CQED 
schemes in which the àtoms are required to be individually addressable p3 . Each atom interacts individually with 
the làser field and the cavity mode, just as for the single atom case. We make the simplifying assumption that the 
àtoms do not interact directly with one another. However, they do experience a second order interaction via the 
cavity mode. The cavity coupling g and the classical pump Rabi frequency r are taken to be idcntical for each of the 
N àtoms, consistent with their indistinguishability p5| . Thus the ./V-atom/cavity interaction Hamiltonian is given, 
for the closed system, by 

N 

where Hi describes the atom/field coupling, Eq. ([ï]), for the ith atom. From now on, we shall use this in the operator 
form 

Hi = -A|1),(1U + r(t) (|l)i(0| 4 + |0)i(l|0 + g (a|l)<<2| 4 + a t |2) i <l| 4 ) , (4) 

where the operators a and af are the annihilation and creation operators for the quantized cavity mode. Since this 
Hamiltonian is symmetric under the exchange of any two àtoms, a symmetric initial state remains symmetric as it 
evolves. We may thus reduce our effective Hilbert state to consider only the states that are completely symmetric with 
respect to atòmic interchange. We may thus use a number representation of the atòmic state, namely |no, ni, r%2, l), 
where n, gives the number of àtoms in state \i) (i = 0,1,2), and l the number of photons in the cavity. We allow 
arbitrary vàlues for no, n\, n 2 and l, and employ the Bose creation and annihilation operators bi and b\ for the atòmic 
states We can then rewrite the interaction Hamiltonian for the closed system as 

H = -Ab\h + r(t) (b\b + ò òi) + g (b\b 2 a + a^&i) . (5) 
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This many-body iV-atom/cavity Hamiltonian conserves the total number of àtoms, represented by the operator 
T = b^bo + b\bi + b\b 2 , as well as the difference between the number of àtoms in state |2) and the number of cavity 
photons, represented by the operator D = b\b 2 — a^a. When cavity decay is added to this description, the operator D 
gains the significance of referring to the number of photons that have escaped the cavity. Since the operators T and 
D commute, we define subspaces e(iV, k) composed of the eigenstates with simultaneous eigenvalues N — «o + Hi + n 2 
and k = n 2 — l of the operators T and D, respectively. This is summarized schematically in Figure ^. 

We find that each manifold e(N, k) contains a null-valued eigenstate given explicitly by 

where Zk is a normalization constant. This state is the analog of the null-valued dark eigenstate for the single- 



atom/cavity system 10 . It contains no population in the |1) infernal state and is thus an iV-atom dark state immune 
to spontaneous decay from any atom. Eq. (||) represents a continuous family of dark states that are transformed into 
each other by the time dependence of r(t). For nearly all vàlues A ^ 0, we can show that the null-valued eigenspace 
of H is non-degenerate, i.e., each of these dark eigenstates possesses a unique energy. We will discuss exceptions 



later (Section 1VC). One can also show that acting on Eq. (Q) with the cavity annihilation operator a, produces 
the corresponding dark state having one less photon, i.e., any such dark eigenstate \tp^) decays to by cavity 

emission. Thus cavity decay does not take the system out of the family of dark states. Conversely, the only way to 
directly couple dark states in different manifolds, is cithcr to spontancously lose a photon or to add a particle. 

This representation suggests that adiabatic evolution might be used for iV-photon generation in an analogous 
manner to that proposed for single-photon generation in Refs. [To). We illustrate this here with the manifold 
corresponding to k — 0. Initially the pump lascr is off, r = 0, and the corresponding initial state has N àtoms in 
the |0) state and no photons in the cavity, (0)) = \N, 0,0, 0). One can then imagine slowly ramping the value of 
the pump làser until the pump lascr coupling is much larger than the atom coupling to the cavity, r 3> g. At this 
final timc, tf, one finds that approximately all N àtoms are in the state |2) and there are N photons in the cavity, 
lV>f(í/))~|0,0,iV,AO. 

This procedure will generate N photons in a closed cavity. We will show that this procedure can be coupled to 
cavity decay to produce an iV-photon state. In practice, realization of this ideal sequence requires that two key issues 
be adequately addressed. First, the need for adiabatic evolution through the family of dark states via couplings 
to the excited atòmic level |1) places constraints on how we vary the strength of the pump pulse, r, based on the 
energy spectrum of the A^-atom/cavity system. Second, spontaneous emission of photons from the cavity will provide 
perturbations to the adiabatic evolution that may be non-negligible. 

We examine these issues in detail in Section IV. Before this, we first analyze the energy spectrum of the closed 
iV-atom/cavity system in the remainder of this Section. This will allow us to establish the critical parameters limiting 
the adiabatic evolution, that are required in order to estimate the errors due to non-adiabaticity and spontaneous 



emission within the open system approach employed in Section [V. While the dark state of each e(iV, k) manifold has 
a succinct description, Eq. (||), the general eigenstates for the JV-atom/cavity system are quite complicated. However, 
in the limit of both large and small r(í), the eigenstates are found to have familiar forms that render them susceptible 
to analytic investigation. 

We examine first the strong pump (large r) limit. Here the eigenstates can be interpreted in terms of the familiar 
angular momentum states. In the limit where r is large relative to the other parameters (g, A), we can neglect the 
terms in Eq. (|^) that are proportional to the cavity coupling parameter g. We then make a transformation from the 
atòmic modes bo and b\ to generalized angular momentum operators, J, using the Schwinger representation. pi], 
This gives J z = \{b\b 1 — bj ) b ),J + = b\b , J_ = b^bi, and J t — h(b\bi + %bo), where J t denotes the total angular 
momentum, J 2 = Jt(Jt + 1)- In this limit we find that the system Hamiltonian becomes 

H = -A(J t + J z ) + 2r(J x ) = -AJ + íl(ij ■ J), (7) 



where Çl = \/4r 2 + A 2 , r\ x = 2r/Q, rj y — 0, and r/ z = —A/Cl. The corresponding energy levels are now identical to 
those of the generalized angular momentum J. Therefore the eigenstates are simply given by \j, m^, Ti2, l) j where j 
is the eigenvalue of total angular momentum, m n is the angular momentum projection along the axis 77, n 2 is the 
number of àtoms in atòmic state |2), and / is the number of cavity photons as before. The cavity coupling term, 

g{b\b· 2 a + ^blbX can now be considered to act perturbatively on Eq. (^), to mix states differing by j — ±|, and 

to change the value of the cavity photon number by unity. Relating these states in the large r limit to our invariants 
T and D, we find that the e(N, k) manifold is composed of states \j, m v , n, l), where < j < (N — k)/2,n = N — 2j, 
and l = N — k — 2j (see Figure ||). Our dark state at large r is |0, 0, N, (N — k))j. In the number representation this 



5 



is simply \nçj — 0, n± = 0, n 2 = N, l = (N — k)) (Eq.(^)). We note that the loss of a photon only reduces the photon 
number l, and conserves j, m rn and n 2 . The eigenstates in the large r limit are shown schematically in Figure (||). 

We now consider the weak pump (small r) limit. Here we find that the eigenstates may also be interpreted in terms 
of a known set of states, but these are now the less well known eigenstates of the Tavis-Cummings Hamiltonian [CT . 
We proceed in this limit by starting from the the system Hamiltonian, Eq. (Q), at r(t) = 0. This is simply 

H = -Ab\b 1 +g(b\b 2 a + a^blb 1 ) . (8) 



We again make a transformation to a Schwinger representation, but this time we choose the transformation to be 
made between modes b 2 and b\ . The generalized angular momentum operators that are created from these two modes 
will be denoted here by F, i.e. F z = \{b\b\ — b\b 2 ), F + = b\b 2 , F_ = b\b\, and F t — ^{b\b\ + b\b 2 ), where F t denotes 
the corresponding total angular momentum, F 2 = F t (F t + 1). Eq. (||) then becomes 

H = -A(F t +F z )+g(F + a + F_aï) ) (9) 

which is recognized to be the off-resonant Tavis-Cummings Hamiltonian ]27| |. Note that this Hamiltonian conserves 
the generalized angular momentum F t . It also conserves sum of the number of photons, l, and the angular momentum 
in the z direction, f z . In the small r limit we can then investigate the effect of finite r using a perturbative analysis. 
This perturbative analysis has two consequences for the energètics. First, the perturbation due to r, which is of the 
form b\bo + òj&i, will only couple states whose total F value differs by 1/2. Second, as a result of this, the resulting 
change in energy of the eigenstates of Eq. (||) is only second order in r. The eigenstates of the Tavis-Cummings 
Hamiltonian are not trivial, but we note that in the limit of large A they are approximately eigenstates of F z with a 
number of photons in the cavity given by l = f — f z — k where / is the total angular momentum and k is an eigenvalue 
of D. The states of the Tavis-Cummings model are described in detail in ]27j and qualitatively in Appendix A. Figure 
|] provides a schematic of these states in the small r limit. 

One uscful advantage of these Schwinger angular momentum representations of the atòmic states for our analysis is 
that in both of these limits of large and small r, i.e., whether for a fixed / or for fixed j, the corresponding eigenvalues 
f z and j z provide a measure of the population in the excited state. This population is given by b\b\ — F t + F z = J t + J z . 



This property will be used in Section IV C to make estimates of the population in the scattering state i.e., in the 



atòmic state |1), that is susceptible to spontaneous emission, and hence of the errors due to spontaneous decay. 



IV. OPEN SYSTEM APPROACH 



The above description of the iV-atom/cavity system as a closed quantum system is clearly incomplete, since a 
proper assessment of the operation of an ./V-photon generator requires the consideration of this CQED device as an 
open quantum system. We must take into account the two channels by which the iV-atom-cavity system interacts 
with its environment. These are: i) the possibility of spontaneous decay from àtoms in the excited state |1) to optical 
modes outside the cavity, determined by the spontaneous decay rate 7, and ii) the coupling of cavity photons to 
electromagnètic modes outside the optical cavity, characterized by the cavity decay rate k |26| . The latter provides 
the required coupling to transfer an ./V-photon state from the cavity mode to an external mode. Cavity decay thus 
plays two different roles here. First, in allowing emission of the ./V-photon state, and, second, affecting the dynamics 
inside the cavity as discussed below. We describe the iV-atom-cavity system as an open system within a quantum 
wavefunction formulation p8| . To characterize its action as an N- photon generator, we evaluate the cavity flux 
2n(a^a). 

In this Section, we use the quantum wave function formulation in a perturbative regime to derive analytic estimates 
for the error rates of its action as an V-photon generator. In Section |v| we then make numerical simulations of the 
full open quantum system that allow us to ascertain the extent of validity of these perturbative error estimates. 

In this open quantum system analysis, we will show that the structure of the closed quantum system, namely 
the identification of manifolds of states e(N,k), each of which contains a dark state that is immune to spontaneous 
emission and that conneets adiabatically to the initial state, still plays a critical role. When r(t) is varied in the 
open quantum system, the N- atom/cavity system evolves primarily within the family of accessible dark states \tp^} 
{k G {0, ...N}) from which no spontaneous emission oceurs, just as in the closed quantum system. The system can 
fail to produce the desired iV-photon output state only if photons are actually lost to spontaneous emission. Thus a 
crucial part of assessing the failure rate of the ./V-photon generator is to quantify the extent of "non-darkness," i.e., 
the probability that the system will evolve towards a bright state from which spontaneous emission may indeed oceur. 
Non-darkness can arise from two factors: non-adiabatic evolution when the rate of change of r(t) is too fast, and the 
conditional dynamics resulting from cavity decay which can couple the dark state \^) to a non-dark state. 
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In the quantum trajectory approach, J^, ^9|, |3(J the dynamics of the ./V-atom-cavity system are characterized by 
dctcrministic non-unitary evolution, interspersed with random 'jumps' determined by photon losses from the cavity 
and by spontaneous emission (scattering) from the excited state. The non-unitary evolution is given by the conditional 
Hamiltonian 

H cond (t) = H{t) - inaïa - i^b\bx. (10) 

The imaginary terms —ina^a and — ijb\b\ describe the back-action on the quantum system from that accumulates 
between instances of cavity decay and spontaneous emission jumps, respectively. The corresponding quantum jump 
operators are given by na and 76. For a full discussion of the quantum jump approach see Ref. p8[ . 

For the numerical calculations prescnted in Section ^ we use the full conditional Hamiltonian, as required in the 
quantum trajectory formulation.^, ^9], |3^| In order to obtain analytic error estimates here, we proceed first by 
assuming that the adiabatic errors are small. Therefore, at all times the state of the system is now regarded as 
deviating only slightly from the dark state. In this situation we may neglect the spontaneous emission term, —ijb\bi, 
treating it implicitly as a higher order perturbation than the cavity decay term. The cavity decay term k is treated 
in first order perturbation theory. Before describing the details of this analysis, we note again that the cavity decay, 
characterized by the jump operator na, connects dark states in the manifold e(N, k) to dark states in a lower manifold 
e(N,k + l), i.e., 



(See Section III above.) Thus neither of the quantum jump operators na or 7Ò1 will lead to errors (i.e. to non-darkness) 
in the operation of our CQED device, as long as the system is maintained within the dark states \^f^). Consequently, 
in order to quantify the failure rate of the ./V-photon generator under these conditions, one need only consider the 
evolution of the system under the conditional Hamiltonian H conc i- We will now make a detailed analysis of the errors, 
starting with an estimate of the extent of non- darkness introduced by cavity decay, then estimating the adiabatic 
error in following the dark state as r(t) is varied, and finally estimating the spontaneous emission flux rate due to 
population in the |1) internal state that is introduced by a combination of cavity decay and non-adiabaticity. 



A. Effect of Cavity Decay on Dark States 

We first analyze the errors due to cavity dynamics. For this purpose, we treat the ina^a term in Eq. (|ï^) as a 
perturbation to the closed-system Hamiltonian, Eq. (||). In particular, under the realistic scenario when cavity decay 
is weaker than the cavity coupling (k < g), we find that the cavity decay term causes the dark states \ip k ) of the 
closed-system Hamiltonian to be modified to the states \ip^} of the conditional Hamiltonian, according to 



w)>= w)) mw) ! — -• ( 12 ) 



Here we have modified our notation to define the states \ip l k (t)) as the instantaneous eigenvalues of the closed-system 
Hamiltonian H conc i with corresponding energies huji. Using this expression, we find the degree of "non-darkness" due 
to cavity decay, e£ av = 1 — |(^ |ï/>°) | 2 , to be equal to 

e ca V = K 2 y m{t)Wam\\ (13) 

From this expression we can generate an upper bound on the degree of non-darkness, e£ av . We first note that 

Y,\(m)Wa\^)\ 2 = ^Kv^WhMtòl'-K^WhMví)! 2 (14) 

= Var^a),^^), (15) 

where Va.r(a' í a)^o^ is the variance in the cavity photon number in the dark state \tp^(t)}. Our upper bound is 
then obtained by replacing the Bohr freqüències u>i in Eq. (|ï^) with the minimal Bohr frequency, and replacing the 
time-varying variance Var(a^a)^o( t ) with its maximum. This results in 

„ max Var(a) a),j,a lt \ 
£ cav < K 2 , ^) (16) 

mm \uij\ 
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The maximum variance of a^a is bounded by N — k (see Appendix B) while the minimal Bohr frequency depends on 
A, and will be discussed detail in Section [V C . Our expression for the bound on the extent of non-darkness due to 
cavity decay thus becomes 



cav < K 2 (N-k) 
mm \u>i\ 



(17) 



B. Adiabatic Errors 



We now examine errors that arise due to non-adiabatic evolution resulting from a non-zero derivative, f(t). Using 
the Standard treatment |3í| we estimate the population leakage from the adiabatic state at a time t, ej^ = | (Q(t)\ïp®) | 2 , 
to first order for a given N and k as 

Here |^(í)) = f t \ipl{t)). We then apply the Standard upper bound for adiabatic error |53|, given by the square of 
the maximum angular velocity of the state, divided by the square of the minimal Bohr frequency, 

£fe " mhTRF ' (19) 

One can bound the maximum angular velocity, max(-0^(í)|-0^(í)) , to be smaller than max(f/ff) 2 (see Appendix 
C). This results in the following upper bound: 



„H (N - k) max(r) 2 , , 

* (k + l)g 2 mm \uj t \ 2 y ' 

Eqs. @ and provide an upper bound for the "non-darkness" in the limit that our perturbative approach 
is appropriate. When the minimal energy separation, min|wj| is independent of N, both equations suggest that the 
"non-darkness" scales at worst linearly with N, in the worst case. Furthermore, we note that the maximum "non- 
darkness" decreases with increasing k. We now proceed to estimate the effect of spontaneous emission, or equivalently, 
the rate of spontaneous emitted flux. 



C. Spontaneous emission flux rate 



The spontaneous emission flux rate is equal to the product of the population in the spontaneously emitting state 
and twice the decay rate of that state, i.e., 2^y(b\bi). Unfortunately, we are unable to analytically calculate the average 
population in the excited atòmic state (b\bi). The simplest way to estimate an upper bound on the spontaneous flux 
rate is then to use the maximum possible value of (b\bi) = N — k, and the maximum probability of the system being 

in a non-dark state, et = e£ av + e^. We can thus obtain an upper bound on the spontaneous emission flux rate as 
2 7 e fc (AT-fc). 

Naturally, by detuning from resonance, we expect to minimize the population that will leak to population in the 
excited state as is the case in both off-resonant Rayleigh and Raman scattering. Hence, we expect the spontaneous 
emission flux rate to change as a function of the detuning A. The value of A aífects the spontaneous error rate in two 
ways. First, it controls the value of the minimal Bohr frequency, and can thereby either reduce or enhance the first 
order population leakage from the dark state. Second, we know by analogy with the three level system (Section || 
and Refs. || |Ï0|) that the "non-darkness" should contain less excited state character as the absolute magnitude of A 
increases. Our scheme works independent of the sign of A. However, for concreteness, A is assumed to be negative 
(red detuning) throughout the rest of the paper. 

Referring to the discussion in Section III, for small vàlues of A, ie.,|A| < r,g, there exist states with energy only 
~ A away from the dark state. In both the small and large r limits, these states are the m — states (Section |ÏTj| ). 
Hence, for |A| < g, the minimal Bohr frequency is equal to A. Furthermore, the eigenstates are cffcctivcly eigenstates 
of J x or F x and therefore have average value J z = 0. Therefore, whatever states the dark state falls into due to either 
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adiabatic error or cavity decay, these states will have large population in the excited atòmic state |1). So in the small 
|A| limit our uppcr bound described above is reasonable. Consequently, we use the "non-darkness" factors calculated 
in Eqs. ( |Ï7| ) and ( p(ï| ) times 2j(N — k). Hence, we find that the error rate for small vàlues of detuning |A| is always 
less than 

/ (jV-fc) 2 max(r) 2 , 2 {N-k) 2 \ 
7 7 ^ (fc + l) 5 2 |A| 2 + |A| 2 )■ { ) 



In contrast, for large vàlues of |A|, the states are effectively eigenstates of J z or F z (Section III). For these states, 
the smallest valued eigenstate is the m — — J or m = —F state, and the corresponding perturbative energy scales in 
second order as E\ = — (fe + l)(.g 2 + r 2 )/A. However, this smallest energy eigenstate has a population in the excited 
state that scales at worse as (fe + l)g 2 /A 2 . We believe that it is a reasonable assumption to multiply this smallest 
eigenstate estimate of the excited state population by our calculated "non-darkness" factor e, in order to arrive at a 
better estimate of the excited atòmic state population. This procedure yields the followng expression for large vàlues 
of A, namely an error rate that is independent of A: 

( {N~k)m^(fY n 2 {N-k) \ 



V (fc + l)V {k + l)g 2 

Careful analysis shows that in the large |A| limit, the state of energy Ei = —(fe + l)(.g 2 + r 2 )/A is not always the 
closest state to the dark state. Indeed states can be identified for particular settings of A and r which are degenerate 
with the null-valued dark state. However, these states are shown to be characterized by higher / or j, and therefore 
in the perturbative limit they do not couple directly with the dark state. Numerical simulations suggest that the 
situation is similar for states which cross the null state in the non-perturbative limit. 



D. Small number of photons in the cavity 



Having obtained error estimates due to cavity decay, non-adiabaticity, and spontaneous emission, we now consider 
the dependence on the cavity photon number l. In the limit that one very slowly ramps up r, the average number of 
photons in the cavity is always much smaller than 1. We assume that our dark state is simply \(N — fc)0fc0) + r)\(N — 

k - l)0(fe + 1)1). For this state Var^a) = rf - -q 4 < rf =< ala >, where rj = ^frr ©)■ In thc limit of 
small ?7, the variance of photons in the cavity equals the number of photons in the cavity. Since the ramp is very 
slow, we expect that the error rate will be due to the cavity dynamics and, therefore, proportional to the variance, 
Eq. (|Ï6|). The rate of photon output is proportional to the number of photons in the cavity. Hence, in the case of 
a small number of photons l in the cavity, we expect at worse for the error to scale proportional to the number of 
photons output. This implies a total error that scales at worse proportional to the number of àtoms N, when A > g. 

One might ask how can we increase r to be of equal magnitude with g, and still maintain a small number of photons 
in the cavity. The reason is that as r increases, there is an increasing chance of a photon being emitted. If one varies 
r slowly enough, the chance of that occurring before r becomes too large is quite high. At this point, we then change 
manifolds e(N, k). Examining the equation for the dark state (Eq. (||)) and looking at only the first two terms, one 
sees that changing manifolds is equivalent to reducing the effective value of r. Compared to the k — manifold, the 

effective r' of the fe manifold is equal to yj ~j^j^j r · As N increases, for small vàlues of fe one then needs smaller and 

smaller vàlues of r for the approximation to hold. For example, when fe = 0, r\ = \fNrjg. However, as the value of k 
increases, r can eventually reach a value comparable to the product Ng and still imply a small number of photons in 
the cavity. 

To summarize, in the limit of a small number of photons l in the cavity and |A| > g, we expect the ratio of 
spontaneous emission flux to cavity flux to be either a constant, or a decreasing function with N. Explicitly, we 

expect spontaneous flux to be smaller than ^2£l JVar(a^a)dt < Jrj 2 dt (see Eq. @ and Eq. ([Ü])). Thc 

ratio of the spontaneous flux, J rj 2 dt, to the cavity flux, 2k J rj 2 dt, is then ^£ which is larger than the fractional 
spontaneous loss, the spontaneous flux divided by the expected number of pho tons out, N. Therefore, for current 
state of the art optical cavity technology for which k = 0.1<? and 7 = 0.05g [[34|, we expect that only 1 of every 200 
photons will be lost to spontaneous emission. 
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V. QUANTUM TRAJECTORY SIMULATIONS 



The quantum trajectory or quantum jump approach allows onc to calculate the properties of open quantum systems 
by averaging over individual quantum trajectòries, p^, [jo| The bàsic elements of this approach were described 
briefly in Section |y|. Technical details of the algorithm are given in Ref. j28| . The calculations described below average 
over 5000 trajectòries per simulation. The simulations were run until the cavity flux and spontaneous emission flux 
were less than 10~ 6 . This results in 2000-5000 steps of length dt = 0.1/ g, depending on the specific parameters of the 
system. 

The bàsic quantity we calculate is the cavity flux, 2«;(aJa), which describes the number of photons emitted from 
the cavity. For these simulations, we have chosen as our figure of mèrit the fractional spontaneous loss. The fractional 
spontaneous loss is defined as the number of photons lost to spontaneous emission, N s , divided by the expected 
number of photons out, N. We calculate N s as the product of 27 and the integral of (b\bi). Due to the statistical 
error resulting from the finite number of trajectòries, we have estimatcd a numerical uncertainty of ±3% in N s /N. 



A. Production of iV-photon state 



Using a simple driving field which increases linearly with time, and realistic cavity parameters based on the cavity 
of Kimble j34j, (f = g/30,K = g/10, 7 = g/20, A = — 2g), we find that the deterministic production of TV-photons 
within a single pulse with small loss is indeed possible. Figure shows the output pulse as a function of time for 
N < 5. We see the expected linear increase in the area of the output pulse, i.e, the cavity flux 2K(a^a), with N. We 
then calculated the loss of photons into the spontaneous emitted modes, shown in Figure ^. When the simulation 
was run in the limit of a small number of photons in the cav ity an d with A > g, we expected that the total photon 
loss will grow at worst linearly, as explained above in Section IV D. Therefore, the fractional spontaneous loss should 
either be constant or decrease with increasing N. This expectation is confirmed by simulations, see Figure ^ and 
Figur e |^. Furthermore, the error is smaller than our expected bound of 0.5% fractional spontaneous loss, see Section 

iyq . 

Although, we have made calculations with only a small number of N vàlues here, it appears that our analytic 
estimates coupled with this numèric evidence for N < 5 suggest that the deterministic production of large N Fock 
states when A ~ g is indeed possible. The numerical results summarized in Figure 6 show that for A > g we have 
minimal fractional spontaneous losses of approximately 0.3actually decrease with increasing N. Since our analytical 
results suggest an upper bound of 0.5expect to be able to produce a Fock state containing on the order of 100 photons 
before the total integrated loss due to spontaneous flux equals a single photon. 



B. Minimizing the Spontaneous Loss 

In the previous subsection, we used current experimental cavity parameters ]34j ] and, for simplicity, a linear driving 
pulse. We now show how one can minimize the spontaneous loss by tuning various cavity parameters and modifying 
the driving field. 



1. Pulse Shape 



For linear ramping with a small rate of change, one finds that the fractional error actually decreases with increasing 
N (Figure §). This is in line with our analytical bounds whic h sugg ested that at worse the fractional error should be 
constant for small occupation of photons in the cavity (Section IV D| ). Figure [j] shows that as the ramp speed increases 
beyond the adiabatic regime, the errors increase rapidly. One can further reduce the error by using a more sophisticated 
time-dependent driving field. For example, the application of a Gaussian pulse, r(t) = gexp — (í — to) 2 / (2r 2 ), reduces 
the fractional spontaneous loss by a factor of 2, relative to the fractional spontaneous loss resulting from a linear ramp, 
as illustrated in Figure |[ One expects that appropriate pulse shaping could further lead to an order of magnitude 
increase in the fidelity. However, the minimal fractional spontaneous loss obtained by pulse shaping is limited by 
errors due to cavity dynamics. In Figure ||[ one sees that this limit is approached by ramp speeds of 5/IOO. 
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2. Effect of Detuning 

The ability to detune from resonance is one of the bàsic tools of atòmic physics . Th erefore, it is important to 



determine how the fidelity of our operation scales with the detuning, A. In Section IV C, our analytical predictions 
suggest that as |A| is increased from A = 0, the initially large fractional spontaneous loss should decrease, eventually 
reaching a constant non-zero value for |A| > g. For larger |A|, our numerical calculations confirm this prediction for 
N = 2 — 5 àtoms, Figure |(| 

In Figure ^, one sees that for detunings, |A|, smaller than the cavity coupling constant, g, the error decreases 
as the detuning increases. One then sees a relatively flat region, followed by an increase in error as |A| increases. 
This rise is outside of the predictive ability of our analytical model and represents a breakdown in the first order 
perturbative methods whcn the energy separation between the dark state and nearest neighbor state becomes smaller 
than the ramping speed. The main conclusion from our numerical simulations is that a wide region of detunings are 
nevertheless possible. In terms of efhciency, (number of photons out per time/fractional spontaneous loss), our data 
suggest that a detuning of |A| ~ (1 — 10)g would be optimal. 



3. The good cavity limit minimizes errors. 



Fo r large vàlues of |A|, one can reduce the error by increasing the value of g. As predicted analytically in Section 



IV C, the total error will reduce as g increases. Furthermore, if one varies g and r proportionally, one can numerically 
observe the 1 / g 2 dependence of fractional loss predicted by Eq. (|23) . This is shown in Figure |l 



VI. CONCLUSIONS 



We have presented a scheme for the deterministic production of V-photon Fock states from coupling of N three- 
level àtoms in a high-finesse optical cavity to an external field. The method relies on adiabatic passage through a 
continuous family of dark states that is controlled by the ramping up of the external field. We have shown that 
this procedure can reliably output pulses of V-photons from the cavity. We have made a detailed analysis of the 
errors involved in the V-photon state production. These result when there is a breakdown of adiabaticity and from 
spontaneous emission. This error analysis yields analytical bounds on the errors which are well reproduced by the 
results of numerical simulations for up to V=5 photons. Our estimates for physically realistic cavity and atom/field 
parameters indicate that this scheme may be used reliably to generate states with up to V=100 photons. 

One way in which such a deterministic Fock state generator may be constructed in the near future is using a 
combination of an ultra-cold atòmic source and a high-finesse optical resonator (with parameters such as those used 
for the calculations in Section |v|) such as those used in the works of Kimble and colleagues p3] . Modern techniques 
of làser and evaporative cooling in magnètic traps can readily produce cold and tightly confined atòmic gases which 
can be transported in vacuum from a production region to the confines of an optical resonator. Using ultra-cold 
alkali àtoms, the states |0) and |2) can be chosen to be two ground state hypcrfinc lcvcls which are both magnetically 
trapped, such as the \F — 1, tuf = — 1) and \F = 2, mp = 1) hyperfine states of 87 Rb. These may be connected by a 
Atof = 2 Raman transition using excited levels on the Dl or D2 transitions - while the strength of this transition is 
suppressed at large detunings |A| from the atòmic resonance, our work shows that only moderate detunings (several 
g) are necessary for high-fidelity operation. The initial V-atom state may be generated by following a weak microwave 
excitation with atom-number-sensitive selection, or, alternately, by real-time observation of the atòmic number in a 
given hyperfine level. This would set proper initial conditions for V-photon generation. 

Several additional comments are in order. Throughout this work, we have assumed that the àtoms are identically 
coupled to the cavity and pump ficlds. Similarly, we assume that these àtoms are indistinguishable in their spontaneous 
emission to optical modes outside the cavity. It would be desirable for the N-photon generator to operate similarly if 
these restrictions are eased. Toward this end, it is still possible to define a family of dark states for àtoms which are 
not identically coupled to the cavity and pump fields. As such, our analytical approach to estimating the photon losses 
to spontaneous emission could be similarly extended to such cases. We believe this would provide similar performance 
to the idealized case we have provided herein, i.e. the deterministic production of Fock states containing as many 
as V = 100 photons should be possible. However, one would expect that while different spatial arrangements of V 
àtoms in a cavity would all produce V-photon Fock states, the specific optical mode occupied by those V-photons 
would differ. Applications requiring many orthogonal, simultaneous pulses of V-photons (such as Heisenberg-limited 
interferometry) would be constrained by these differences. 

We envisage a number of applications for deterministically-produced V-photon states. As alluded to in the above 
discussion, two orthogonal V-photon states can be combined using linear beam-splitters to create an interferometer 
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which is sensitive to differential phase shifts between the two arms of the interferometer which are on the order oíl/N 
(the Heisenberg limit), rather than the typical 1/yN sensitivity (the Standard quantum limit) obtained with classical 



light pulses 1 36 . Two orthogonal iV-photon states can also be used to create a highly entangled state using only mea- 
surement and linear òptics p7| | . Such highly entangled photon states can be used to perform precision measurements 
pj[ |39| , ESI. The CQED device described in this work could thus be used for demonstrations of this interferometric 
method. While the low photon numbers deterministically available from this CQED device (given current state-of-the- 
art cavity parameters) would not yield the precision available from the use of much more intense classical light sources, 
there may be applications rcquiring high precision at low light levels which are enabled by this Fock-state generator, 
e.g. coherent control fH]] . Other applications in the field of quantum information, communciation, cryptography, 
and computation are also possible, and we intend to investigate such potential applications in future work. Finally, a 
deterministic N-photon source would find both bàsic and applied use for the absolute calibration of optical detectors, 
particularly those designed to be sensitive to múltiple photons (as opposed to conventional avalanche photodiode 
devices). Similarly, the A-photon generator can be used as the light source for novel multiphoton spectroscopy, e.g. 
REMPI Q, although the optical frequency range which can be produced by a given atom-based system would be 
quitc limited. 
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VIII. APPENDIX A 



In this appendix, we present a qualitative picture of the eigenstates of our system in the small r limit, i.e., a 
description of the eigenstates of the non-resonant Tavis-Cummings Hamiltonian, Eq. (Q). First, note that the cre- 
ation/annihilation operators associated with the cavity were not present in Eq. (B), the Tavis- Cummings Hamiltonian 
would become again simply a Hamiltonian describing angular momentum (F) about some axis (z). In fact, for large 
enough A, one sees that is indeed the case, and we do merely have eigenstates of angular momentum about the z 
axis. In order to get a qualitative picture of the Tavis-Cummings eigenstates for general vàlues of the cavity coupling 
g and finite vàlues of l, it is useful to treat the cavity creation/annihilation operators as operators that act to enhance 
the "effective magnètic field" in the "i" direction. For example, when one has a large number of photons, l, in the 
cavity a Standard approximation is to replace the cavity photon creation/annihilation operators with \fl. Applying 
this transformation to Eq. (JoJ) , one finds that the Hamiltonian becomes simply H — A(F t + F z ) + g\/lF x , i.e., the 
cavity coupling has been transformed into an effective magnètic field in the x direction. 

The eigenstates of Eq. (^|) are constructed as follows. First, one starts with atòmic states that are the eigenstates 
of F z , namely, |/, f z ) = \ni) = f + f Z) n-j = f — f z . We know that the Hamiltonian conserves D, the difference in the 
number of àtoms in state |2) and the number of photons, in the cavity. In this representation, the conserved quantity 
is —f + fz + l = —k. We can then append a photon state that refleets this conservation symmetry to the atòmic states 
of modes |1) and |2) identified via F z . This results in the atom-photon states \f, f z , 1í)f — \ f, fz) f\1 — ~k + f — f z ). 
The only missing component is now an atòmic state of mode |0). We then append the appropriate number state of 
the atòmic mode |0) such that N is conserved, according to N = 2/ + nç,. This then specifies the remaining index no. 
We have thereby arrived at a set of states \f, f z , k, N) = |n = N — 2f, n\ = f + f z , ni = f — f z , l = — k + f — f z ), 
that may now be làtex npb.used as a basis for expansion of the exact eigenstates of Eq. (^|). Since the Hamiltonian 
Eq. (|^) conserves /, k, and n, the exact eigenstates can be formally written as |/, m, fc,n)j? = c m f z \f,f z ,k,N). 
Naturally, determining the coefficients c m k is not trivial p7| . For our purposes here, it is not necessary to find explicit 
solutions for the coefficients. We require only the main features of the energy spectrum, which are available from the 
identification of the above basis and are described below. 

It is important to note that in order for the photon number l to be positive, we must have / — f z > k. For fixed / 
and k vàlues, this can lead to a truncation of the possible accessible states. One finds that for positive k, the angular 
momentum projection, m, can have only 2(/ — k/2) + 1 vàlues, whereas for negative k, m can take all 2/ + 1 vàlues. 
This implies that for k < 0, the eigenstate degeneracy for given vàlues of /, N, and k, equals the degeneracy one 
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would expect for a total angular momentum value of / (i.e., 2/ + 1). In contrast, for vàlues k > 0, the eigenstate 
degeneracy equals the degeneracy expected for a total angular momentum value oi f — k/2, i.e., 2(/ — k/2) + 1. Note 
that in our scheme we always assume that there are initially no photons in the cavity, i.e., k > 0. 

In summary, we see that one think of the states \f,m,k,N)p as being cffcctively "angular momentum" statcs 
possessing a total angular momentum F t equal to / — k/2 and angular momentum projection of m about the "magnètic 
axis" . Although this analogy is not perfect because of the spread of actual eigenstates over this basis, it does contain 
the following important feature. This is, that for A larger than g\Jl ma x — gV%J — k, where l max is the maximum 
number of photons in the cavity, the eigenstates are eigenstates of F z , with eigenvalue m — k/2. So the large A limit 
is simplified. In contrast, for small and medium size vàlues of A, the eigenstates are superpositions over a wide range 
of f z states. 



IX. APPENDIX B 



In this Appendix, we derive an upper bound on the variance of cavity photon number for the dark state. In order 
to reach this bound, we need to first outline other properties of the dark state. We begin with a few definitions. The 
dark state in the E(N, k) manifold is given by, 



N-k 



l**> = i- y ( X) \(N -k-l) (l + k) l), (23) 

' ' Z k £jj y/(N-k-l)\{l + k)\V. 



=0 

where x — r/g and Z k is the normalization constant 

d 

dx 



Let be the state given by -f-|*fc)- Thus, 



1 / 1 N ~ k K-tY- 1 7' \ 

where Y k is the normalization constant. 

We further define the state |xfe) given by the normalized state H\^> k )- Therefore, 

N-k-l ,_ w 

\Xk) = T7r Y \(N -k-í-l) (l + k) l), (25) 

IA ' W k y/(N-k-l-l)\(l + k)\ir 

where W k is the normalization constant. 

We now present some properties relating these normalization constants to each other and their derivatives. So let 
(l)k be the average number of photons in the cavity for the darkstate in the k manifold. In other words, 

(l)k = (* fe |a t a|* fe ) 

N-k 



1 J ^ Ix 21 

zïL· {N -k-iy.(i + k)m 

Property 1 



3 = (l)k (27) 
7 k 



Proof: By definition, Z k is the normalization constant of the dark state \^k}- Hence, 

N-k 

^ 2 = E 



N ; x 21 



(N — k — l)\(l + k)W. 



.: 2Z' k Z k = 



2lx 2l ~ 1 



.: xZ' k Z k = J2 



(N - k - 1)1(1 + fc)M 
Ix 21 



(N - k - l)\(l + k)\l\ 
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Dividing both sides by Z 2 gives us, 

N—h 



xZ' k 1 Ix 2 



Z 2 X 



Z k Zl j^{N-k-l)\[l + k)\l\ 
= (l) k [FromEq. f 



Property 2 



Zk+i x 
x- 



Z k 

Proof: Z k +i is the normalizaton constant for the state j^fc+i). Hence by definition 



= (l)k (28) 



(N- (k + 1) - + k + 1)U\ 

(iZfc+i) - 2^ (jv-fc-i - /)!(/ + fc + i)m 



i=0 



, 7 ,2 ^ (i + l)- 2(i+1) 



i=0 
iV— fe) 



(JV-fe-(í + l))!(í + l + fe)!(í + l)! 



(.xZ fc+ i) 2 = ^ 



2/ 



^ (7V-/fc-Z)!G + fc ) m 
Dividing both sides by Zf , we get the result, 



(xz k+l v = i y ) (/)x 2 ' 

V ^ / ^ 2 (iv-fc-0!(i + fe)m 



From Eq. fl2q), we prové the required result. 
Property 3 

If the system is in the dark state \^ k ) : then 



Var(a)a) = {xY k f (29) 



Proof: It is trivial to see that Property 3 is true when x = 0. Let x =/= 0. Consider the state \& k ) given by eq. U4 
Since i/Owe can extend the summation to l — 0. Thus, 



1 / 1 N ~ k U-tV" 1 7' 



i jV-fc - , . / \Z-1 

r fe Z fc ^ V 2*/ V(^- fc -0!0 + 



[Using Eq. |23| 

1 ï=í r-rV- 1 

à >/(jv- jfc-z)!(i + fc)m' 

[Using Property 1] 



E^-Wfc) 7^ = = i == l(iV-fe-Z) (í + fc) o 
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Yfe is the normalization constant of |$&). Thus, 

Y " = ^ S ^ _ {N-k-l)\(l + k)\l\ 

1 x 21 

: - {xYk)2 = zlg^-ffl fe ) 2 (Jy - fc -0!(f + fc)m 
.-.(^n) 2 - (/ 2 ) fe -((0fe) 2 

It easy to see that this is merely Var(a' r a). 
Property 4 

(l 2 )k = (l)k((l)k+i + l) (30) 

Proof: 

{P)k = ZÏ § (i\r-fc-Z)!(/ + À:)!Z! 
Jv-fc - ,2; 



3? ^ 

j iV-Jfc- 

= z 2 



Ix- 



Z 2 Z^ {N -k-l)\{l + k)\(l-l)\ 
(l + l)x 2 V+V 



z l U {N-k-i-m+i+k)\i\ 

2 2 N-(k+l) 



(AT- (fe + l)-Z)!(Z + fe + l)!Z! 

x Z fe+1 J_ + ij^ 



z fe j=Z (iv-(fe + i)-0!(í + fe + i)m 

= -2^(/ + l)fc+l 

Using Property 2 we get the required result. 
Property 5 

(l)k > (l)k+i (31) 

Proof: We need to show that 

(l)k - (Ofc+i > (32) 

which is equivalent to showing that 

Z 2 k Z 2 k+1 (l) k - Z 2 Z 2 k+1 (l) k+1 > (33) 

since Zj is positive for all j. 

N-k , 9; N-k—1 2m 

Ix 1 o \ - mx 



z ï z k+i( l )k - z%z% +1 (i) k+1 - X! (m_u_ JWi -+- tv/ï ~ ^ 51 



^ (N — k — + k)\l\ fc ^ o (iV-A-l-m)!(m + ife + l)!m! 

JV-fc-l 2m ^-fc 7„,2i 

Ex ^ — ^ tx 
T/V — i- — 1 — mWm 4- h 4- lli'm! 



(JV - k - 1 - m)!(m + fe + l)!m! ^ {N-k- + k)\l\ 
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N—k 9/ N — k — 1 9™ 

X x -> mx 



E (N - h- l\Hl 4- h\W E 



(N - k - + k)W. ^ (N - k - 1 - m)!(m + k + ï)\m\ 
(N — k — m)x 2 



r 2(N-k) i"-]?- 1 (N-k- mk 2 " 



E 



N\(N - 1 =0 (N-k-l- m)!(m + k + l)!m! 



We note that 



Thus, 



(N - k - 1 - to)! (to + fe + l)\m\(N - k - + k)U\ 

- 2{N - k) Ty^ 1 {N-k- m)x^ \ 

Nl(N-k)\\^ Q [N - k - 1 - m)\(m + k + l)\m\ J ~ ' 1 ' 

^ +1 <0fc-^+i<0fc+i> E E (]v _ k - 1 - m)!(m + k + l)\m\(N - k - /)!(/ + jp ' (35) 

In order to evaluate the sum in Eq. (^||), we choose two integers between and N — k—1, a and b. We then evaluate 
the sum of the two terms corresponding to l = a, m = b, and l = b, m = a. We see that 

(0 - b)x 2(a+b ï (b - a)x 2 ^ a+ ^ 



(N-k-l- b)\(b + k + l){b\(N -k- a)\(a + k)\a\ (N-k-l- a)\(a + k+ ï)\al(N -k- b)\(b + k)\b\ 

x^2(a + b)) (, (N-k-b N-k-. 

' (a — b) 1 



(N — k — b)l(b + k)\b\(N — k — a)\(a + k)\a\ \ '\ b + k a + k 

x l 2(a+b)) ^ n „;„,„!„ „„„J + „ J t OTTr ,;„ ^q+V^t- ( n _ M f N-k-b _ N-k-a 



Since {N _ k _ mb+ky ^ {N li_ a) , (a+k) , a] > 0, we simply need to determine whether (a - b) [-^ ^fct J > 0. If 

a > b, a — b is positive and ( N ^ k b — ^tï+fc~ a ) ^ s P os itive, so the product is therfore also positive. If b > a, a — b is 

negative and (^ N b ~^ k b — N a+k~ a ^j IS a ^ so negative, so therefore the product is still positive . Hence, for all a and b 

2 ( a+fc ) (, / N-k-b N-k-a\\ , , 

-TTT >0, (36) 



(JV- fc- 6)!(& + fe)!6!(iV- fc-a)!(a + fc)!a! V V 6 + fc a + fc 

which implies that 

N-k-l N-k-l W 2(m+/) 
Ü èo - ^ - 1 - + A; + 1)!ttz!(A^ - fc - Z)!(Z + A;)!Z! ~ ' 1 ' 



By Eq. Q35J) and Eq. (|37|), we have 

Z 2 k Z 2 k+1 (l) k - Z 2 k ZÍ +1 (l) k+1 > (38) 

and thus 

(0* > (Ofc+i- (39) 

Property 6 

If the system is in the dark state |\^fc), 

Var(a)a) k < (N - k). (40) 

Proof: 

Var(a)a) k = (l 2 ) k - ((l) k ) 2 

= (l) k ((l) k+l + í)-((l) k ) 2 
[Using Property 4] 

< (l)k 

[Using Property 5] 

< N-k. 
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X. APPENDIX C 

We determine here an upper bound on the angular velocity of the dark state, maxa 2 = max($>(t)| $?(*)). We 
first note that using the chain rule the m&xa 2 is equivalent to determining the maximum value of x 2 Y 2 where Y^ 

is the normalization constant defined in equation |24|). For convenience, we will define c; = = x ^ and 

H <— v ' (N -k-l)\(l+k)\l\ 

i — \ l 

c[ = ^ k 7) ! (z+fc)!» (^ nrou S nou ^ this Appendix, we will use /' to represent the derivative of / with respect to x, 

and / to represent the derivative f with respect to í) . 

Since the maximum of the product of two functions is always less than or equal to the product of the maximums 
of each function, the first step is to simply take the maximum of x 2 . This yields 



, 2 

maxr 



maxi 2 = [—y-j ■ ( 41 ) 

The second step is to determine the maximum value of Y k 2 . Using the notation from Appendix B, we first note that 
(*fc|$fc) = y-^kl^k) = °' since I**) is normalized. From Eq. Q), this implies that 

N-k 

(* k \Y,4\(N-k-l) (l + k) l) = Z' k . (42) 



1=1 



Hence, we can write the following expression for Y k 



2 




y2 = 4s[YM) 2 -W\- <^> 



Since Y k 2 must be positive and both X)/=i ( c í) 2 an d (Z'k) 2 are positive, we can write the following inequality 



r2 ^ S/ = 1_( C Í) 2 

zl 



Y 2 < " . (44) 



J k 

Note that Y k 2 is an even function of x, so the value of Y k at x — must be a local minimum or maximum. When one 
takes the derivative of Y/c with respect to x, one finds only a single zero at the origin. The identity, x 2 Y k = var(cva), 
implies that Y k must go to zero as x increases in order to maintain a finite variance. Therefore, the maximum value 
of Yk occurs at the origin. Explicitly calculating the limit at the origin yields 

N — h 

maxY 2 = ±— -. (45) 
fe k + 1 

Consequently, we have an expression for the maximum angular velocity 

2 (N — k) maxf 2 

maxai. = — . (46) 

(k + l)g 
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FIG. 1: Pictorial representation of the Hamiltonian for a single three-level atom in a single mode cavity. The atom is driven by 
an external driving field of frequency ui r and coupling strength, r. The atom is also coupled to a cavity mode of frequency ui c 
and coupling strength g. Both the cavity and the driving fields are detuned from the atòmic transition resonance by a common 
frequency A. The atom/cavity states are denoted here as \i, l) where i = 0, 1, 2 are the three atòmic levels and l is the number 
of cavity photons. 
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FIG. 2: Pictorial representation of the Hamiltonian for 5 àtoms in the e(5,3) and e(5,2) manifolds written in the number basis 
representation and assuming red detuning, A < 0. The e(N, k) manifold is composed of the eigenstates with simultaneous 
eigenvalues TV = no + n\ + ni and k — n^ — l where rn is the number of àtoms in atòmic state i and l is the number of photons 
in the cavity. The transition strength between levels is proportional to either the driving field, r, or the cavity coupling, g. 
The energy of the states in the absence of all couplings is given by -niA. For red detuning, A < 0, the states with a higher 
occupation of the exited atòmic state, ni , will have a greater energy. The dark state is the superposition of states with ni = 
described in Eq. ([]). A transition from the manifold e(5,3) to e(5,2) occurs when a photon is emitted from the cavity. This 
transition preserves the dark state (see text). 
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FIG. 3: Pictorial representation of the Hamiltonian for 5 àtoms in the e(5,3) manifold using the Schwinger representation for 
atòmic modes |0) and |1) and assuming red detuning, A < 0. This representation is appropriate when the external driving 
field, r, is much larger than the cavity coupling, g. The eigenstates to íirst order are eigenstates of angular momentum about 
an axis defined by the effective magnètic fields B z — — A and B x = r. The cavity coupling, g, acts as a perturbation, mixing 
states differing by a total Schwinger angular momentum of 1. 
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FIG. 4: Pictorial representation of the Hamiltonian for 5 àtoms in the e(5,3) manifold in the Tavis-Cummings basis for atòmic 
modes |1) and |2) and assuming red detuning, A < 0. This representation is appropriate when the cavity coupling, g, is much 
larger than the external driving field, r. The eigenstates of the Tavis-Cummings Hamiltonian are complicated (see text and 
Ref. [^|) but preserve total Schwinger angular momentum / obtained from modes |1) and |2). The coupling to the external 
field, r, acts as a perturbation, mixing states differing by a total Schwinger angular momentum of 1. 




FIG. 5: The cavity flux, 2k{o í a), for N = 2 — 5 àtoms plotted versus dimensionless time, T = gt. These simulations were done 
using a linear ramp r(T) = RT where R = p/100. The other cavity parameters were 7 = g/20, k — g/10, and A = — 2g. The 
integrated cavity flux, N s , increases linearly with iV as expected. 
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FIG. 6: Fractional spontaneous loss for N = 2 — 5 as a function of A. A linear ramp was used with R — g/100, 7 = g/20 and 
k — g/10. As expected one sees that in the large A limit that the fractional spontaneous loss scales less than N (see text). Our 
analytical bounds in the large |A| limit suggest that the fractional spontaneous loss should be less than = 5 x 10~ 3 (dotted 
line). The discrepancy at large |A| represents parameters for which our perturbative approach is invàlid. For the small |A| 
limit the increasing fractional spontaneous loss is consistent with leakage to states that have a higher occupation of the excited 
state. 
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FIG. 7: Fractional spontaneous loss as a function of ramp rate, R, and detuning A (7 = g/20, and k = g/10). As expected, 
larger ramp rates lead to an increased loss when the system can no longer adiabatically follow the dark state. Note how the 
dependence on A varies with R. In the adiabatic limit, small R, large A yields the lowest íractional spontaneous loss. However, 
when one ramps the system quickly, i.e., at large R, the dependence on A dependence is reversed and small A yields the lowest 
loss. 



24 



3- Exp 



e- Lin 




2 1 2 3 4 5 

N 

n 

FIG. 8: Fractional spontaneous loss for linear and Gaussian ramping (A = —2g,j = g/20, and k = g/10). For the linear ramp 
(circles), R = g/30. For the Gaussian ramp (diamonds), r(t) = gexp(t — ío) 2 /(2r 2 ), r = 50/g. The Gaussian width, r, was 
chosen such that the process of emitting a single photon would occur in the same time as the linear ramp. One expects that a 
more sophisticated pulse could result in an order of magnitude reduction of the fractional spontaneous loss. 
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FIG. 9: Fractional spontaneous loss as a function of cavity coupling, g, and detuning, A (R = g/100, 7 = k/2). The vàlues oí 
A in the legend are in units oí 10«. One sees that at all vàlues of A there is a significant decrease in fractional spontaneous 
loss for increasing g. In the large |A| limit, one fi nds th at the decrease in spontaneous loss with g has a í/g 2 dependence, as 
predicted by our pertubative analysis (see Section IV D). 



